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"^v^". We study the cosmic microwave background (CMB) temperature fluctuations non-gaussianity due 

I ■ to the vector mode perturbations (Alfven waves) supported by a stochastic cosmological magnetic 

| field. We present detailed derivations of the statistical properties, two and three point correlation 

£\j , functions of the vorticity perturbations and corresponding CMB temperature fluctuations. 
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I. INTRODUCTION 

In the framework of the standard cosmological scenario the cosmic microwave background (CMB) temperature 
1 fluctuations are gaussian and are fully determined by the CMB temperature fluctuation two-point correlation functions 
while higher order odd correlation functions (for example three-point correlations) are identically zero (for a review 
Q_[ on CMB fluctuations and possible non-gaussianity see Ref. [lj and references therein) . This is a consequence of the 
inflationary scenario that predicts the gaussian initial perturbations, and even at the level of perturbations there is 
no violation of rotational symmetry 1 . Several observations indicate that the CMB temperature map could be slightly 
non-gaussian |3| and thus to adequately describe the CMB fluctuations one must study higher order correlations 
functions. Furthermore, some modifications of standard inflationary models lead to a slightly non-gaussian CMB map 
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A common way to characterize the CMB temperature fluctuations non-gaussianity is to introduce thefNL parameter, 
which in fact determines the relation between the two-point and three-point correlation functions [5(. The current 
CO ! CMB data limits / NL = 32 ± 21 Q . PLANCK mission will be able to give us with stronger limits (to improve current 
Tj" ■ limit by an order or a few). If the nearest future CMB measurements, for example PLANCK data will confirm that 
|/nl| is not significantly less then one, the standard cosmological scenario must be revised substantially. There are 
several different ways for such a revision. For example the inflation could be driven by multiple fields, or the Universe 
isotropy has been violated at very early epochs [lj . Recent studies [f| 0) address the magnetic field induced density 
perturbations as a source of the CMB temperature fluctuations non-gaussianity. The physical meaning of this effect 
is as follows: the temperature anisotropics caused by the magnetic field are proportional to the magnetic field energy 
density parameter, AT/T cx B 2 /p CI Q, where p cr is the critical density today and B is the comoving value of an 
■ effective magnetic field. Accounting that the square of the magnetic field B 2 is a non- linear form Q, the corresponding 
AT/T fluctuations must be non-gaussian. The limits for the scale invariant magnetic field amplitude from the CMB 
rjj . non-gaussianity test is of order of 1CU 9 Gauss @. 

In this paper we investigate the CMB non-gaussianity due to the vector mode of perturbations induced by a 
stochastic magnetic magnetic field, (see Refs. |10l - [l2l ] for details of the vector magnetized mode). If the magnetic field 
presence is a reason for the CMB non-gaussianity, this magnetic field must satisfy several conditions: (i) the magnetic 
field must be generated in the early Universe, prior to recombination; There are different mechanisms to generate 
magnetic fields in the early Universe, such as inflation, phase transitions, see for reviews [IH; (ii) the correlation 
length of the cosmological magnetic field must satisfy the requirement of causality [T3 |. and thus to have a field 
correlated over the horizon scale or even larger this field should be generated during the inflation and have a scale 
invariant spectrum |15l | ; (iii) the amplitude of this magnetic field should be small enough to preserve the isotropy of 
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1 The generation of the vector mode, which involves a preferred direction, is not excluded during the inflation, but the exponential 
expansion washes out the vector (vorticity) perturbations if no supporting external source is present 
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the background Friedman-Lemaitre-Robertson- Walker (FLRW) metric (so the energy density of the magnetic field 
should be the first order of perturbations) , be below the upper bounds (few nGauss) given by observations [H] . On 
the other hand the amplitude of the magnetic field should be large enough to leave observational traces on CMB. 
Recently it was argued that non-observation of blazars in TeV range by Fermi mission indicates the presence of large 
scale correlated intergalactic magnetic field with a lower bound of order of 10~ 16 Gauss fl71| . The existence of the 
lower bound of order of 10~ 16 Gauss magnetic field favors the magnetic field of cosmological origin [18], and thus the 
magnetic field amplitude at 1 Mpc is squeezed between 10~ 9 and 10 -16 Gauss. 

Vector and tensor modes of magnetized perturbations are much more complicated then the scalar one. The first 
paper to address the CMB bispectrum induced by the vector and tensor mode of perturbations has been Ref. [l9| 
where the analytical expressions were derived. A natural extension of that approach, namely numerical estimation of 
the vector and tensor modes induced non-gaussianity were presented in Refs. [l9|-|22[. In our study we give detailed 
derivations of the two and three point correlation functions of the CMB temperature fluctuations induced by the 
magnetized vector mode. In this sense at this stage this article has a methodological nature. 

The outline of the rest of the paper is as follows: In Sec. II we define the magnetized vector mode of perturbations 
and compute the Lorentz force two- (Sec. IIA) and three- (Sec. IIB) point correlation functions. We explicitly 
discuss in details the difference of the vorticity perturbations two-point correlations functions, and show that the 
non-gaussianity of the vector field (vorticity) can be seen already from the two-point correlation function. In Sec. Ill 
we address the CMB temperature fluctuations induced by the magnetized vector mode perturbations. We present 
analytical expressions for the two- (Sec. IIIA) and three- (Sec. IIIB) point correlations of the CMB temperature 
anisotropics. We briefly discuss our results and conclude in Sec. IV. Useful mathematical formulae and details of 
computations are given in Appendix. 



II. MAGNETIZED PERTURBATIONS VECTOR MODE 



To study the dynamics of linear magnetic vector perturbations about a spatially-flat FLRW homogeneous cosmolog- 
ical spacetime background (described by the metric tensor g^ — a 2, q^ v , where 77 M „ = diag(— 1, 1, 1, 1) is the Minkowski 
metric tensor and a(rj) the scale factor) we follow the standard procedure and decompose the metric tensor into a 
spatially homogeneous background part (g^ v ) and a perturbation part, g^ = g^ u + 5g^ v , where /i, v £ (0, 1, 2, 3) are 
spacetime indices. Vector perturbations 5g^ u can be described by two three-dimensional divergence-free vector fields 
A and H where 

$goi = %o = a 2 Ai, Sg i:j = a 2 (H id + H j:i ). (1) 

Here a comma denotes the usual spatial derivative, i,j £ (1,2,3) are spatial indices, and A and H vanish at spatial 
infinity. Studying the behavior of these variables under infinitesimal general coordinate transformations one finds that 
V = A — H is gauge-invariant (the overdot represents a derivative with respect to conformal time). V is a vector 
perturbation of the extrinsic curvature [23| . Exploiting the gauge freedom we choose H to be constant in time. Then 
the vector metric perturbation may be described in terms of two divergenceless three-dimensional gauge-invariant 
vector fields, the vector potential V and a vector representing the transverse peculiar velocity of the plasma, the 
vorticity ft = v — V, where v is the spatial part of the four- velocity perturbation of a stationary fluid element [ll|. 2 
As we noted in the Introduction, in the absence of a source the vector perturbation V decays with time (this follows 
from the equation for V, V + 2(a/a)V = 0) and therefore can be ignored. 

The residual ionization of the primordial plasma is large enough to ensure that magnetic field lines are frozen into 
the plasma. Neglecting fluid back-reaction onto the magnetic field, the spatial and temporal dependence of the field 
separates, B(£,x) = B(x)/a 2 [251 ] . Since the fluid velocity is small the displacement current in Ampere's law may 
be neglected; this implies the current J is determined by the magnetic field via J = V x B/(47r). Accounting for 
a frozen- in magnetic field lines the induction law takes the form B = V x (v x B). As a result the baryon Euler 
equation for v has the Lorentz force L(x) = — B(x) x [V x B(x)] /(in) as a source term. The photons are neutral so 
the photon Euler equation does not have a Lorentz force source term. The Euler equations for photons and baryons 
are [H [ij] 

Ji 7 + f(v 7 - v fc ) = 0, (2) 



2 Given the general coordinate transformation properties of the velocity field v, two gauge-invariant quantities can be constructed, the 
shear s = v — H and the vorticity CI = v — A [231 . In the gauge H = (i.e., V = A) we get O = v — V l'24l . 
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a „ r 



L( v \x) 



n h + -Sl b - p(v 7 - v&) = — - — r, (3) 
a R a(pb+Pb) 

where the subscripts 7 and b refer to the photon and baryon fluids, and p and p are energy density and pressure. 
Here f = n e ara is the differential optical depth, n e is the free electron density, <jt is the Thomson cross section, 
R = (pi, +Pb)/(p-y +P7) — 3 j Ob/4 j o 7 is the momentum density ratio between baryons and photons, and is the 
transverse vector (divergenceless) part of the Lorentz force. In the tight-coupling limit v 7 ~ v;,, so we introduce the 
photon-baryon fluid divergenceless vorticity f2 (= 0~ = Clt) that satisfies 

(i + R)ii + R-n= 4 L(y)(x) , . (4) 

a a 4 (p 7 +p 7 ) 

The average Lorentz force (L(x)) = — (B x [V x B]}/(4tt) vanishes, while the r.m.s. Lorentz force (L(x) • L(x)) 1 / 2 is 
non-zero and acts as a source in the vector perturbation equation. 

To proceed one needs to obtain an expression for the Lorentz force in terms of the magnetic field characteristics. 
We assume that the magnetic field is Gaussian and satisfy 3 , 

(B 4 (k)B;(k')) = (2Tr) 3 S(k-k')P lJ (k)M(k) for k < k D , (6) 

and vanishes for k > kr>. Here a star denotes complex conjugation, and <5(k — k') is the usual 3-dimensional Dirac 
delta function, and ko is the magnetic field damping scale defined through the Alfven waves damping [l(| We 
approximate power spectrum M{k) by simple power laws. 

Following Ref. [ll| it can be shown that the Fourier transform of the Lorentz force l\ V \]s.) = fcIL(k) is related to 
the Fourier transform of spacial part of magnetic field energy momentum tensor 7^ . (k) 



T$ \k) = ^tAs J d 3 q[B 4 (q)^(k - q) - ^S l3 B m (q)B m (k - q)] , (7) 



as 



Hi(k) - i(F y -(k)fc m + ^(k)^)^ (k) , (8) 

where Py(k) = 5y — kikj is the transverse plane projector with unit wavenumber components ki — ki/k. 

Next we shall solve the Eq. Q. It can be solved in two different regimes, for length scales larger and smaller then 
comoving Silk scale As- These two solutions are (llj : 

* For the scales A > As, (k < ks) 

(i + i?)(p 70 +p 10 ) 

where p l0 and p 7o are photon energy and pressure today. 

* For smaller scales with A < As, (k > kg) 

n(M)= ,, r rc (k) ^ — T> (10) 

where L 7 is the comoving photon mean-free path. 

Eqs. © and ([TUf define the vorticity through the magnetic source, n(k), Eq. ©. As we can see the vorticity 
perturbations at scales below the Silk damping stay constant, while the vorticity perturbations above the Silk damping 
are increasing linearly with the time. In order to compute the magnetized vector mode induced effects we need to 
compute the vorticity perturbations two- and three-point correlations. 

Below we present our computations. It must be stressed that the form of the magnetic field two-point correlation 
function, Eq. ([6]) presumes the following properties of the field: (i) transverse, divergence free field, V • B = 0, and 
in the fourier space this is insured by the projector Py(k); (ii) isotropy (no preferred direction) insured by S(k — k') 
and axis-symmetry of the field, the two point correlation is symmetric under under i and j indices exchanges; (iii) 
The magnetic field is a gaussianly distributed field. 



For a vector field F we use 

F J (k) = J Ae^FjW, F J (x) = J ^Le-^FjOO, (5) 

when Fourier transforming between real and wavenumber spaces; we assume flat spatial hypcrsurfaccs. 
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A. Lorentz Force Two-Point Correlation Function 

For the two-point correlation function of the II (k) magnetic vector source we have 

C^>,k0^n*(k)n i2 (kO>. 

Using results of Sect. I A 21 one finds for it 

C| l 1 f 2 ) (k,k') = 5(k-k')VlS(k), 

with 



(11) 



(12) 



(8tt) 2 



d 3 qM(|q|)M(|k-q|) 



P iia (k)k b + P Hb (k)k a P l2C {k)k d + Pi 2 d(k)fe c Pac{q)Pbd(k - q) + P ad {q)P bc (k - q 



(13) 



Note that as usual we assume summation over the repeated indices. Using -Pjj(k) projector symmetry properties and 
defining j = k • q, j3 = k • (k — q), and fi = q • (k — q), after long but simple computations we obtain for the Lorentz 
force two-point correlation function 



(LP(k)Lf(k')) = ^g^5(k-k')/rf 3 qM(|q|)A^(|k-q|)(2-/3 2 - 7 2 



(8tt)2 



+ ^S(k - k') / d\M(\d\)M(\k - q|) [2 7 2 /3 2 - 7 2 (1 - /3 2 ) - /3 2 (1 - 7 2 )] 



(87T)' 



5(k - k') / d 3 qAf(|q|)Af(|k - q|) - /? 2 ) + (k - q) 4 (k - q),(l - 7 2 ) 



j(hqj + ^©5(1 - 2/3 2 ) - /3 [fc(k - q)j + fe,-(k - q) 4 J (1 - 2 7 2 ) 
7/3 q t (k - q)j + q 5 (k - q) f } . 



(14) 



Note that the form of the Lorentz force two-point correlation function, Eq. (|14j) . has several symmetries: (i) symmetry 
under i and j index exchange; (ii) the symmetry under q and k — q exchange (i.e. 7 and /3 angles exchange symmetry). 
It must be underlined that the trace of Eq. (fT4"|) is given by 



(iP(k)if ) (k')) = -^5(k-k') J d 3 qAf(|q|)Af(|k-q|)(l + M7 /3-2/? 2 7 2 ) 



(15) 



which totally agrees with the result of Refs. [1, [Til E3- Further simplification of Eq. (|14|) gives following result for 
the Lorentz force two-point correlation function: 



(Lp(k)if »(k')) = J^^(£)J(k-k') | d 3 qA/(|q|)A/(|k-q|)(l + A17 /3-2/3 2 7 2 ) . 



(16) 



The ensemble averaging procedure insures that rotational isotropy is preserved. Note that without ensemble averaging 
in one realization the isotropy of the two point correlation can be violated and only the averaging leads to cancelation 
of anisotropic component and restoration of isotropy. 



B. Lorentz Force Three-Point Correlation Function 

The three-point correlation function of the vector magnetic source is given 4 by 

Cl^ki^ks) = (n il (k 1 )n i2 (k 2 )n i3 (k 3 )) . 



(17) 



4 It is appropriate to define the three point correlation function as an average of three n^'s without the complex conjugation. 
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Using results of Sec. I A 21 one finds 

(ki , ka, k 3 ) = *(k! + k 2 + k 3 )^7 3 (ki , k 2 ) , (18) 

with 

VffiO^ki.) = l^qMdqDMdki-qDMdka + ql) 

^tiji(ki)Pj ia3 (kx - q)(ki + k 2 ) a 3Pi M - 2 (k2)Pj„- 3 (k2 + q;/', , ik, + k 2 )fei sl P slS2 (q)fe 2s2 • (19) 
The Lorentz force three point correlation function is given as 

(4 V) (k 1 )i^ ) (k 2 )Li 3 y) (k 3 )) = JLj(k 1+ k 2 + k 3 ) | ^qMdqDMdki-qDMdka + ql) 

Pm (ki)P 2 ^ 2 (k 2 )P 3i3 (k 1 + k 2 )P, lS3 (k 1 - q)(ki + k a ).,P iai ,(ka + q) 

[(kx X q)(k 2 X q)] . (20) 

We see that the r.h.s. of Eq. ([2D]) is symmetric under exchange ki and k 2 . It is obvious that this symmetry is 
reflected in the CMB temperature three point correlation function, see below. 

III. TEMPERATURE FLUCTUATIONS FROM THE MAGNETIZED VECTOR MODE 

Vector perturbations induce CMB temperature anisotropics via the Doppler and integrated Sachs- Wolfe effects 

mm, 



(x ,n,r,o) = -v-n|^ cc + / dr,V ■ n, (21) 

J ??dec 



T 



where n is the unit vector in the light propagation direction, ?7 dec is the conformal time at decoupling. Xo denotes 
the observer position, x doc = xo + n(r;o — 77 dec ). Due to the spherical symmetry the temperature fluctuations are 
decomposed using the spherical harmonics as, 



AT v 



(x ,n,7? ) = ai m {x ,r] ) ■ Y lm (n) . (22) 

1=0 m—-l 

Accounting for the spherical harmonics property, see Chapter 5 of Ref. [29| . we obtain az m (xo,77o) = 
J dil n l^* n (n)AT/T(xo, n, ?7o). The decaying nature of the vector potential V implies that most of its contribu- 
tion toward the integrated Sachs- Wolfe term comes from near 77 dec . Neglecting a possible dipole contribution due to 
v today, from Eq. ([21]) we obtain, 

AT 

— (n) ~ v(?7 dcc ) • n - V"(r? dec ) • n = 0(77 dcc ) • n . (23) 

Here we placed observer at the origin Xo = and we skip 770 denoting AT/T(n) = AT/T(xo = 0, 11,77 = 770) and 
ai m = a-imfa-a — 0, ?7o)- Since x dcc = 11(770 — ?7dec) Fourier transform of Eq. (|23p results in 

AT 

— (k, n) = n(k, Vdcc ) ■ n)e < < k - n > A " , (24) 

where r2(k, 77 dcc ) is the Fourier amplitude of vorticity perturbations at r\ — ?7 doc , wave vector k = /ck labels the 
resulting Fourier mode after transforming from the coordinate representation x to the momentum representation by 
using e ikx , and A?/ = 7/0 — ?7 dec ~ 770 is the conformal time from decoupling until today. 
For the multipole coefficient Fourier transform we obtain, 



ah 



,(k) = J dn(tt(k,77 dcc ) • n)e^Y^(h) , (25) 
where we use a; m (k) ee a; m (k, 77 =770). 
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A. Two-point correlation function 

The two-point correlation of the temperature fluctuations is given by 

AT AT 

C(n 1 ,n 2 ) = ( — (m) — (n 2 )) = X, IJ ( a hm 1 ahm 3 )Y^ mi (n 1 )Y l2m2 (n 2 ) , (26) 

l x mx l 3 m 3 

here ^2 tm = X)z*^o Y^ m =-l- ^ n t ne case wnen the spatial isotropy and rotational invariance is preserved the multipoles 
with different I and m do not correlate and (a* imi ai 2 m 2 ) = Ci 1 Si 1 i 2 S mi m 2 - It is obvious that C(ni, n 2 ) is the function 
of the angle between ni and n 2 vectors, and we can easy get 

07 1 1 

C(n 1 ,n 2 )=J2^ r C l P l (n 1 -n 2 ). (27) 

Irn 

The sound waves (density perturbations) sourced by the magnetic field keep the rotational invariance unchanged 
and thus the multipole correlation matrix has a diagonal form, and the two-point correlation function depends on 
one angle, ni ■ n 2 . The same time the temperature fluctuations are non-gaussian @, 0| due to the non-linearity of 
the magnetic field energy density with respect to the magnetic field. Note, that the magnetic field which has a 
non-gaussian energy momentum tensor Q can be itself field with Gaussian distribution. 

As soon as the vector mode is present in the Universe there is an additional direction inserted through the vorticity 
field r2(?7doc). Before ensemble averaging procedure the two-point correlation is not rotationally invariant and the 
cross correlation function between the multipoles has off-diagonal components [13, [H| ■ The measurement of these 
off-diagonal terms might serve as a tool to constraint the primordial homogeneous magnetic field. 

We define the multipole two point correlation function as usual: 

f <i^k d^k 

D hl 2 m2 = ( a hmx a hm 2 ) =J 7^I72^p"( a ? : imi( kl ) a i2m2( k 2)) , (28) 
and the simple calculation leads to 

^^dn^e'^"^^^ .(29) 

We also decompose the vorticity perturbation plane wave over the vector spherical harmonics as, r2(k)e lk n,?0 = 
i m ^-im^im ( n )' w ith ^ = —1,0, 1 and use V • tt = condition, leading to fi(k) • Yj" 1 ^ (k) = 0, see also Appendix 
of Ref. [2S| . We use also the definition n • ^^(n) = Yj~^*(n). Then we have 



KC 2 = I |^S^ n ^m>l)*W , (30) 

fin t 2 r 2 



= W- 1 V^TT)^(^(k) • Y^k)) . (si) 



with 

The integration over n in Eq. (|30p gives us Si ltl , Si 2t2 , S miri , S m2r2 , and the sums are eliminated and we get 

D nr im2 = (-1)' 1 -' 3 i h+l '- 2 y/h(h+ l)h(h + 1) f d 3 k 3h(k lVo )j h (k 2 Tlo) 

hh (2ir 2 ) 2 J 1 2 fcife?? 2 



Y^r(ki). Y^ (k 2 ) . (QK^nM). (32) 

Next we have to consider separately case of large scale approximation Eq. Q and case of small scale approximation 
Eq. ([ID)). Using Eq. (JIBJ we obtain: 

(2^)^2(1 + i?dcc)(p 70 +Pj )\ \ Vo J J 

(1 + M7/3 - 2 7 2 /? 2 ) (y<£>*(£) • Y<+« (k)) (33) 
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for the scales larger then the Silk damping scale. For the scale smaller than the Silk damping scale we have 



- " (2 7 r)3[2(i 7 /5)(p 70+P70 )P yu — *»viH,,-v.- HU (%))2 

(1 + M7 /3 - 2 7 2 /3 2 ) (y£>*(£) • Y^ 2 (k)) . (34) 



Recall that we assume that the magnetic field power spectrum is given by a simple power law, M(|q|) oc q n . To 
proceed we have to evaluate the following integrals over angular variables 



q (l + M/?7 - 2 7 2 /3 2 )(fc 2 + g 2 - 2fc 97 )"/ 2 (Y^f(k) • (k)) . (35) 



It is easy to see that integrals over dSlg. and df2 q separate. First we evaluate the integral over dSl q . For a particular 
k, choose the polar axis for the dO q integral in the z direction. Since 7 — cos #3 the integrand does not depend on 
the azimuthal angle <f)q and the integration over cfta simply gives 27r. The integration over dcos# q is simple to be 

evaluated and is given as j\ d 7 (l - 7 2 ) (l - k l q ^^ ) {k 2 + q 2 - 2fcg 7 )" /2 , with q = |q|. Thus 

2™ = 2nS hh S mim2 f d7(l - 7 2 ) f 1 - il 7 S!!lL ) ( fc2 + « 2 - 2fc 97)" /2 ( 36 ) 



fc 2 + q 2 - 2fc<? 7 

and finally only the diagonal cross correlations are present. Note that, the rotational symmetry (absence of off-diagonal 

-jmi' 



terms in /D™ 17 ™ 2 ) is due to the ensemble averaging. 



B. Bispectrum definition and calculation 

The standard approach to study the CMB non-gaussianity consists on the CMB temperature fluctuations inves- 
tigation. Below we present the self-consistent way to describe the bispectrum. As usual the three-point correlation 
function of CMB temperature anisotropy is defined as 

AT AT AT 

£(ni,n 2 ,n 3 ) = ( — (m) — (n 2 ) — (n 3 )) = ^(a hmi ai 2m2 ai 3m3 )Y hmi (n 1 )Yi 2m2 (n2)Yi 3ms (n 3 ) ,i = 1,2,3 , (37) 

li17li 

In order to proceed we need to calculate the following form 

/d 3 ki (i 3 k 2 d^k 3 
(27r) 3 (27r) 3 (27t) 3 ( a 'i m i( kl ) a '2™2( k 2)a; 3 m 3 (k3)) • (38) 

Using expression for a; m given by Eq. (|2B1 we get 

d 3 ki d 3 k 2 d 3 k 



\^/imi ^Jg-ma ^3 777-3 / 



I d ki d k 2 d k 3 i(k 1 .n 1+ k 2 .n 2+ k 3 .n 3 )A 7? 
7 (2^) 3 (2tt) 3 (2tt) 3 1 2 3 

Y l imi ( n 2)^*™ 3 ( n 3)«H 1 «2 l2 ^3 (Oil ( k i . m C c)Sli 2 (k 2 , 7J dec )fi i3 (k 3 , 77dcc)) • (39) 

First we integrate over dxii by using Eq. (117) on p. 227 [2^. Proceeding in the way given in Sec. IIIA we arrive at 

B iM 3 - (2^)3 ViiCi + l)fc(*a + l)ia(l3 + 1) y d kxd k x d k 3 F^T^ 

Y^f (kOlnY^f (k 2 )| i2 Y^*(k 3 )| i3 (k 1)?7dec )Q i2 (k 2 , 77 dcc )^ 3 (k 3 , %cc )) . (40) 



Next we have to consider separately case of large scale approximation Eq. Q and case of small scale approximation 
Eq. (HO}. Using Eqs. ([17 ] [15 1 1X5 ]) . relation Eq. (TA"8j) and the representation for the Dirac delta function Eq. (TA"9|) after 
some computations for the large scale approximation (L > Lg) we obtain 



^ = (2 7 r 3 ) 3 [(l + fl doc )( P70+P70 )] 3 [HTJ J^^^^^M^M^) 
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dSl tl dSl ta dSl ia J dfi q Jq 2 dqM(\q\)M (|k x - q|)Af(|k a + q| 

x 2 dxj tl {kix)j t2 (k 2 x)j t3 (k 3 x) G'i* 



EEE*' tl+t2+t3 

*H*1 *2r2 t 3 T-3 



"1*1*2*3 



YgWxq Y^ 2 (k 2 )x(k 2 + q 

J a L 

Y t ^(k 3 ) X (ki - q)] ft [Y z « 3 *(k 3 ) X (k 2 + q) 
where Glj. 1 ***? is the usual Gaunt integral, 

^nra'ra = / ^Y tm (A)^ 2 , 2 (n)Y t3 r 3 (A) 



Y^^kOxCkr-q) 

(i) *, 



(41) 



(42) 



The indices a, 6, c correspond to the vector components and repeated ones reflect the scalar product of the corre- 
sponding vectors. 

Similarly for the small scale approximation (L < Lg) we obtain 



j^m 1 m 2 m 3 
hl2h 



dk 1 dk 2 dk 3 j h (k 1 r] )ji 2 (fc 2 r?o) ji 3 (k 3 r) ) 



(2 7 r3)3[(L 7 /5)(p 7o+ p 7o )]3 ?70 3 
dQ^dn^dQ^J dftq J g 2 d g M(|q|)M(|ki-q|)M(|k 2 +q|) 

E E E %tl+t2+t3 J x2dx ^ (Mj* ^ G r\ 



•>*1*2*3 

ir 2 r 3 



*nTi t 2 r 2 t 3 r 3 
Y t ( - 1 2 ) (k 2 )xq 



x q 



Y i ( 1 1 i/(k 1 )x(k 1 -q) 



r( l ) 
' lim 2 



(k 2 ) x (k 2 + q) 



Y t ( 3 - 1 3 ) (k3)x(k 1 -q) 



,(1) * 



(k 3 ) x (k 2 +q) 



(43) 



Using the Wigner D functions and proceeding in the way analogous to Sec. IIIA it can be shown that there are 
non-zero cross correlations between non-equal / and m. Again the proper answer assumes accounting for the angular 
dependence in the power spectra M. 



IV. CONCLUDING REMARKS 



We present the study of the magnetized perturbation vector mode induced CMB two- and three-point correlation 
functions in a common framework. We show that already the vorticity two-point correlation functions reflect the 
anisotropy of the considered perturbations before ensemble averaging. One of the implications of our results is that 
CMB bispectrum computation technique presented in Ref. (30j in the case of magnetized perturbations should be 
applied with caution. 

Note that in the present work we have focused on derivation of main analytic results for the CMB two- and three- 
point correlation functions arising from the vector mode supported by the stochastic cosmological magnetic field. We 
plan to present a detailed analysis of obtained equations and phenomenological estimates in a separate publication. 
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Appendix A: Useful Mathematical Formulae 

In this Appendix we list various mathematical results we use in the computations. 

1. Wigner D Functions 

Wigner D functions relate helicity basis vectors e'±i = =F( e e ± ie^)j\f2 and e'o = e r to spherical basis vectors 
e±i = T{ e x ± ie y )/V% and e = e z (sec Eq. (53), p. 11, [29j]) through 

e M = E^(^ e '°) e -' i/, At = -1,0,1. (Al) 

V 

In both the spherical basis and the helicity basis the following relations hold: e u e^ = 5 V/X , e M = (— l) M e_ Al , e M = e*, 

1. Calculation of B 

For calculation of the two-point correlation function we must determine the magnetic field energy momentum two- 
point cross correlation {t^P (k)T^f^ (k')) given by Eq. A6 of 11 1 and Eq. 4.1 of 8]. It is easy to show that the parts 

of the t^P (and rffi proportional to the 6 c d), see Eq. ([7]), do not contribute to the integral of Eq. <j8j) . Then we need 
to compute the following object 



*Wki,k 2 )= J (^3 < g a(qi)^(ki ~ q 1 )fl e (q 2 )^(k 2 - qa)) 



(A2) 



Using Eq. ^ and Wick's theorem, we obtain that the contribution of the two-point correlation of the magnetic field 
energy momentum into the vorticity perturbation is given by 



<5(k-k') 



J d 3 qM(\q\)M (|k - q|) [p ac (q)P bd (k - q) + P ad (q)P bc (k - q)] . (A3) 



(4tt) 2 

For calculation of bispectrum we need to know following object 

*W e /(ki,k 2 ,k 3 ) = j ^^^(^(qi) B "(ki -qi)P c (q 2 )P d (k 2 -q 2 )B e (q 3 )i? / (k 3 - q 3 )) . (A4) 

Assuming that the magnetic field obeys Gaussian statistics we can expand the six-point correlation function using 
Wick's theorem. Doing this we will get seven terms proportional to either 5(ki), 5(k 2 ) or <5(k_3) (which we neglect) 
and eight terms proportional to 5(ki + k 2 + k 3 ) which we keep. So, for B a bcdef one finds (§1. l3~i|: 

Ba6cde/(ki,k 2 ,k 3 ) = <5(k!+k 2 +k 3 ) J d*qM{\q\)M{\\<L X -q\)M{\\<L 2 + q\) 

{P oc (q)n e (ki - q)P*/(k 2 + q) + P, (q)A/(kx - q)P de (k 2 + q) 
+ Pad(q)P 6e (ki - q)P e /(k 2 + q) + P od (q)P 6/ (k 1 - q)P ce (k 2 + q)} 

+ 5(k!+k 2 +k 3 ) J d 3 qM(|q|)M(|ki-q|)Af(|k 3 +q|) 

{P ae (q)P 6c (k 1 - q)P d/ (k 3 + q) + P ae (q)Ad(ki - q)P d /(k 3 + q) 
+ P af (q)P bc (k 1 - q)P de (k 3 + q) + P a/ (q)P 6d (k 1 - q)P ce (k 3 + q)} . (A5) 

Assuming that projector A abcde f acting on this object is symmetric w.r.t. each pair of indexes, i.e. 

j^abcdef j^bacdcf j^abdcef j^abcdfe (A6) 

quantity B a bcdef can be brought to a compact form 
B a 6 cde /(ki,k 2 ,k 3 ) = 8d(ki +k 2 +k 3 ) J d 3 qM(|q|)M(|ki - q|)M(|k 2 + q|)P ac (q)P 6e (k 1 - q)P d/ (k 2 + q) . (A7) 
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3. Useful relations 

Following relation from vector algebra is valid (see e.g. Eq. (31), p. 16 [29(): 

[A x B] ■ [C x D] = (A ■ C) (B • D) - (A ■ D) (B • C) . (A8) 

Using plane wave decomposition into spherical harmonics one obtains following useful representation for the Dirac 
delta function 

J(k! + k 2 + k 3 ) = J 7^e*( kl + k2+k3 ) x = 8 J d^^^^** 1 ^^ 3 ^^!^^^^) 

tiri t2T2 

y tiri (k 1 )y* ri (x)y t2r2 (k 2 )y4 2 (x)y i 3 r3 (k3)Y t * 3r3 (x) . (Ay) 
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